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Abstract. We investigate late time acceleration of the universe in higher dimensional cosmology. The 
content in the universe is assumed to exert pressure which is different in the normal and extra dimensions. 
Cosmologically viable solutions are found to exist for simple forms of the equation of state. The parameters 
of the model are fixed by comparing the predictions with supernovae data. While observations stipulate 
that the matter exerts almost vanishing pressure in the normal dimensions, we assume that, in the extra 
dimensions, the equation of state is of the form P oc p 1-7 . For appropriate choice of parameters, a late time 
acceleration in the universe occurs with qo and z tr being approximately —0.46 and 0.76 respectively. 
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1 Introduction 

By now, several observations have confirmed that our universe is in a phase of accelerated expansion. Fur- 
ther, it is also clear that this acceleration had set in at a relatively recent time in the cosmological calen- 
dar. Starting with observations of Type 1A supernovae[l, 2], this feature of the Universe is now strongly 
suggested by observations of cosmic microwave background (CMB)[3], large-scale structures [4], baryon 
acoustic oscillations[5] and gravitational lensing[6]. Although observational confirmations are piling up, a 
convincing theoretical framework is still lacking. Several interesting mechanisms have been suggested to ex- 
plain this feature of the universe, such as cosmological constant[7], quintessence^ ], modified gravity[9-l 1], 
chaplygin gas[12, 13] and many others. However, these models have their own shortcomings. For exam- 
ple, models with a non-zero cosmological constant need a high degree of fine tuning [8, 14-19] whereas 
potentials required for quintessence models are unnatural in the context of particle physics[20]. Appeal- 
ing to higher dimensional cosmological models is another promising mechanism to explain this mysterious 
phenomenon. This is the line of approach we adopt in this paper. 

There are various possible constructs in the extra dimensional context including (but not limited to) 
brane world models [21-23]. Here, we consider a particular simple model akin to that used in Ref[24]. 
Whereas Ref [24] invoked extra dimensions to solve the horizon problem in early universe, using an anisotropic 
fluid residing in 1 + D\ + D2 dimensions, we adapt the formalism to produce a late time acceleration instead. 

Motivated by observations, we assume that the universe is filled with a uniform density matter. How- 
ever, the pressure exerted by the matter in the normal dimensions is different from that in the compact 
dimensions, while being isotropic within each subspace. As we will argue, observations severely constrain 
the functional dependence of the pressure on the density. Within this constraint, however, a very simple 
form of the equation of state gives an excellent agreement with data. Although pressureless matter would, 
normally, decelerate the expansion of the universe rather than accelerating it, it is the inteiplay with the 
hidden dimensions that provides the impetus for this expansion. 

It should be clarified at this stage that our construct is not a brane-world scenario and that we do 
not attempt to address issues such as the hierarchy problem in the Standard Model of Particle Physics. 
We, rather, make the simplifying assumption that these extra dimensions are compactified to a scale small 
enough to play essentially no direct role at the TeV scale. In this sense, the spirit is closer to more canonical 
scenarios defined in dimensions larger than four (an example could be a generic model derived from String 
Theory). Possible phenomenological manifestations of such models are postponed to future discussions. 

The rest of the paper is constructed as follows. In section 2, the formalism of the model is developed, 
and the equation of state argued for. In the subsequent section, we present the solutions to the ensuing 
evolution equations. In Section 4, we compare the predictions of the model with data and infer the preferred 
values of the parameters. And, finally, we conclude in Section 5. 
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2 Evolution Equations 



We start with a spacetime which has, in addition to one temporal and three normal spatial dimensions, D 
extra spatial dimensions. This 1 + 3 + D dimensional spacetime is described by the line element 

ds 2 = _ dt 2 + a 2 [t) (-^— 2 + r 2 d^j + b 2 (t) (j^p + R^n-l) • (2.1) 

We reserve the super(sub)script '0' for the time dimension. Whereas lower-case Roman indices = 
1,2,3) denote the normal spatial dimensions, upper case Roman indices denote the extra dimensions and 

take the values I, J = 4, 5, , D + 3. Here, D is a parameter which takes integral values and is to be fixed 

by comparing with observations. 

In eq. (2.1), a(t) denotes the scale factor in the normal (3)- dimensions and b(t) represents the scale 
factor in the extra dimensions. Since we consider the entire (1 + 3 + L>)-dimensional universe to be homo- 
geneous, the two scale factors a and b are functions only of the time t. As is well known, the visible universe 
is well described by a vanishing spatial flatness (k± = 0) and we shall assume the situation to be so. For 
reasons of simplicity as well as symmetry with the observed sector, we shall assume = as well. 

For this line-element, the components of the Einstein tensor (with p,,u = 0,1, ... ,D + 3) read 

a b a A 2 tr 

Ct--*-^-^ V, (23, 

a b a b a z 2 b z 

G\ = (1 - D)\ - 3^ + 3(1 - D)±\ + (D - 1)(1 - £■ £ - 3^J V/ . (2.4) 
b a ab 2 b z o z 



The energy-momentum tensor is assumed to be of the form 

T$ = diag(-p, P a , P a , P a , P b ,..., P b ) (2.5) 

where p is the energy density of the fluid and P a (Pf,) is the pressure exerted in normal (extra) dimensions. 
This form of energy momentum implies that there is isotropy within the subspace associated with the nor- 
mal dimensions and also within the orthogonal subspace spanned by the extra dimensions. However, the 
pressures in the two subspaces are different. Note that it is the observed large-scale isotropy of the universe 
that prompts one to consider an isotropic matter distribution. No such restriction applies to the pressure 
exerted in the extra dimensions, and thus, we could as easily have considered more elaborate structure for 
Tj. However, other than adding more freedom to the model, this would not have resulted in any particular 
qualitative improvement to the scenario. Hence, we desist from adopting such a course and adopt eq. (2.5). 
The fact that T„ needs to be divergenceless (T^ M = 0) implies 

±(pa 3 b D ) + P a b D U* + P b a^b D = . (2.6) 
at at at 

In standard (3 -dimensional) cosmology, the constituents of the universe today are dark energy (essentially in 
the form of a cosmological constant), dark matter and baryons. (The radiation energy density has substan- 
tially redshifted and hence, has negligible contribution to the energy density of the universe.) Of these, both 
dark matter and baiyons are very well approximated by pressureless matter. Indeed, the nature of the dark 
matter can be inferred very well as compared to the dark energy. Consequently, we start by assuming that 
the matter does not exert any pressure in the three visible directions and P a = 0. On the other hand, in the 
extra dimensions, it does exert a pressure of the form P b = P b {p)- Note that this implies that it is a single 
fluid that exerts such an anisotropic pressure. A mechanical analogy would be that of gas molecules filling 
a space but constrained to move only along a subspace. We could, of course, have adopted a scenario with 
two fluids, each inhabiting a subspace. However, this would only have increased the degrees of freedom 
in the theory without adding qualitatively to our understanding. Hence we desist from doing this, although 
such a course of action may well be necessitated when one attempts to construct a microscopic theory. 
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The dark energy is, thus, directly 'visible' only in the extra dimensions. The effect (late time accelera- 
tion) in the normal dimensions is through the evolution of extra dimensions. An infinite variety of equation 
of states for this dark energy are possible. For simplicity, we will assume a monomial form, viz. Pf, = w^p 
with 

w 

w b = — , (2.7) 

P 1 

where w and 7 are parameters of the model to be chosen so as to reproduce the observational data and p is 
defined in eq. 3.1. Once again, this choice (reminiscent of a generalized Chaplygin gas [25]) also serves to 
minimize the number of free parameters in the theory. 

The aforementioned energy momentum tensor governs the evolution of the 1+3+D dimensional space- 
time. Let us now consider the Einstein equations, namely 

G$ = nTjf, (2.8) 

with k = 8ttG. The '00' component can be expressed as 

d 2 hb D(D-l)b 2 8nGp 
a z a b o tr 3 

Solving for a fa and bringing it to a form close to the familiar form of the FRW equations, we have 

d 2 8irGp D (2D + 1) b 2 Db Id (D + 2) b 2 32 nGp .„ 1m 

^ = ^T + 6 ¥ T 2bxi 3 V 2 + ^T~ (2 - 10) 

In the absence of the last two terms on the R.H.S., this equation would, understandably, reduce to the 
standard form. In other words, the b/b dependent terms act as an effective dark energy source 1 . 

Since the scale factor, b(t), corresponding to the extra dimensions enters the equations only through 
b/b (= dlnb/dt), what is relevant for the evolution of a(t) is not the absolute value of b but only the ratio 
by which b(t) changes with time. This is because in our model, spatial curvature in the hidden dimension 
&2 = 0. In other words, it is not the size of the hidden world that matters, but its fractional rate of change 
(compression or expansion). It is this rate of compression of the extra dimensions that effectively acts like a 
dynamical dark energy source for the visible universe. Note that, for a non-zero b/b, the scale-factor in our 
world, a(t), evolves non-trivially even in the absence of any matter (p = 0). This is not unexpected, because 
gravity does couple the two subspaces and the contraction (expansion) of one can lead to the expansion 
(contraction) of the other. To be specific, p = leads to a power-law evolution of the two scale-factors (the 
exponents being determined by D) wherein one of them increases with time and the other decreases. 

As is well known, usual (1 + 3)-dimensional deSitter cosmologies admit both expanding and contract- 
ing solutions. We choose one of the solutions, namely, the expanding one, because we observe that the 
universe is expanding and not contracting. The situation here is a little more subtle. As far the '00' com- 
ponent of the Einstein equations goes, there is still a generalized symmetiy of the form (d, b) (— d, —b). 
However, as the form of the other components of shows, this symmetry is not manifest. A consequence 
of this and the preceding discussion is that a pressureless fluid (even in a higher-dimensional world) would 
not admit 'late time acceleration', although uniformly accelerating/decelerating solutions are possible[26]. 

3 Cosmological Solutions 

The evolution of the universe is governed by the Einstein equations along with eq. (2.6 & 2.7). Not all of 
these are independent, though. For example, using the constraint equation (G® = 8irGp) and the continuity 
equation (T^ M = 0) we may eliminate b. Before we do so, it is convenient to rescale the variables in terms 
of dimensionless quantities, namely 

t .. a' a 

f = A = = 

Hn a a Hn 

(3.1) 



P = — 



Pc b bH { 







'it is worth pointing out that, in models of non-minimal coupling, terms that are linear in a/a do appear on the right hand side 
of the FRW equations as in eq. (2.9) 
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where l/p c = 8ttG/(3Hq) and primes denote derivative with respect to r. In terms of these variables, the 
equations of motion now read 

= (D + 2) A" + 3(D + 1) A' 2 + D (1 ~ D ^ B' 2 + D (D - 1) A' B' + 3 D n pw b 



p = -p [3 A' + D (1 + w b ) B'] 
B' = (D- l)" 1 



(3.2) 



3A'± J3D- 1 <j (D + 2) A' 2 + 6 ^— - ^ pfi 



These are two first order differential equations in A' and p with the last of the three being an algebraic 



(3,-3.10, 0.66) 
(4, -2.62, 0.64) 
(6, -2.80, 0.74) 
, -2.50, 0.60) 




1.4 
1.35 

1.3 
1.25 

1.2 
1.15 

1.1 
1.05 
1 



(3,-3.10, 0.66) 
(4, -2.62. 0.64) 
(6,-2.80, 0.74) 
(8, -2.50, 0.60) 



-0.6 -0.4 
I 



Figure 1. Behaviour of the scale factors a(r) [left panel] and 6(r) [right panel] with the rescaled time r. The numbers 
in the parentheses refer to (D, w, 7). 

relation. Note that the Einstein equations can only determine A'(t) and B'(t) and not the scale factors 
themselves, a situation exactly analogous to the (1 + 3)-dimensional case. To solve eq. (3.2) for A'(t) 
and p(r), we require two initial conditions. Since we know the conditions in the present universe relatively 
precisely, we will prescribe the conditions today. In other words, with r = referring to the present epoch, 
we evolve these equations back in time with the following 'initial' conditions 



a 
a 

P_ 

Pc 



A'\ T=0 = 1 



T=0 



T = 



(3.3) 



P\r=0 



We may now numerically solve the two coupled first order differential equations. Two such solutions exist, 
one for each sign in the last of eq. (3.2). We reject here the branch with the '— ' sign as it leads to an 
accelerated expansion for all times rather than a transition from a decelerated phase to an accelerated one. 

The model is characterized by three parameters namely D, 7 and w. In figure 1, we present the solu- 
tions for the two scale factors for some representative values of these parameters. For ease of comparison, 
we have rescaled the solutions 2 so that 



«|r=0 
b\r=0 



A\ T -_ 

BL-_ 








We have a whole class of solutions in which a, the scale factor for our universe starts from 1 at r = 
and decreases monotonically for negative values of r. A word of caution is in place here. Since we have 
neglected radiation completely, the equations are valid only as long as the universe is matter dominated. The 
redshift corresponding to matter-radiation equality is about z eq ~ 2.9 x 10 4 f2 m /i 2 . For the measured values 



This is not to say that the two scale factors are indeed the same in the present epoch, but reflects the fact that the scale factors 
are arbitrary upto a constant. 
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of Q m ,(~ 0.27) and h(~ 0.72), we have z eq ~ 3000. In fact, even before we go as far back as z = z eq , the 
approximation breaks down as the radiation density can no longer be neglected. We have checked, though, 
that the inclusion of the radiation component does not change the evolution drastically for z J> z eq . And, 
since we are primarily interested in the evolution of the universe in relatively recent times, the inclusion of 
radiation does not affect the results in any discernible way. 

As can be easily discerned, the relative evolution in b(r) is small. In the time interval that a(r) has 
increased by nearly a factor of 3000, 6(r) has decreased by ~ 29%. While the opposing signs of the 
evolution was predicted even in a matter-free universe (see discussion in the preceding section), the large 
difference in the magnitude of the evolution is but a consequence of the difference in the pressure exerted 
by matter in the two worlds. 

The epoch of matter-radiation equality, r eq = r(z = z eq ), has a considerable dependence on the 
parameter choice. For a given value of D and 7, a smaller w shifts r eq further into the past [for example, 
with (D, 7) = (6, 0.59), we have r eq « -0.95, -0.97, -0.99 for w = -2.58, -2.6, -2.62] thus increasing 
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Figure 2. Evolution of the deceleration parameter q with the redshift z. The numbers in the parentheses refer to 

(D,w,i). 

the present-day age of the universe. Similarly, for a given value of D, and w, as we decrease 7, once again r eq 
shifts further into the past [(D, w) = (6, -2.8) leads to r e<? = -0.91, -0.88, -0.86 for 7 = 0.7, 0.72, 0.74]. 
Note, though, that the ranges of parameters are restricted (and correlated). An arbitrary set would tend to 
destroy either late time acceleration or shift r eq to unacceptable values. 

As figure 1 shows, the curves for a(r) have a slight upward concavity for r ^ —0.3. This is but a 
reflection of late time acceleration. Prior to this epoch, the universe was in a decelerated phase, as attested 
to by the prominent upward convexity at r ^ —0.5 . This becomes clearer when we plot the deceleration 
parameter q as a function of the redshift (see figure 2). 

4 Observational Constraints 

Having established that the model, for some choice of parameters, does lead to correct late time acceleration, 
we now seek to confront it with other observational data. The most important such data relates to Type la 
supemovae. The very comprehensive Union2 data set [27] lists the distance modulus fi as well as the redshift 
for 557 such supernovae. As the distance modulus fi = 5 logd^ + 25 is nothing but a rephrasing of the 
luminosity di{z), defined through 

z 

d L (z) = (l + z) f^- y (4.1) 


we, then, need to calculate di(z), given our determination of H(z) for a particular choice of parameters. 
We define a x 2 -test through 
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where /x t ^ defines the value expected in our model for a particular choice of parameters, whereas /x & s and a 
are the observational value and the associated root-mean-squared error. We may, now, determine the best-fit 
value of the parameters by minimizing the x 2 - I* 1 table 1, we list such best fit values for some choices of D. 
The results are analogous for other choices. Note that, with an increase in D, both \w\ and 7 decrease. This 
is quite understandable as the ensuing smaller extra-dimensional pressures would now have an enhanced 
effect in the normal world owing to the larger effective coupling between a(t) and b(t). Thus, if we were to 
admit veiy large D values, without any concern for the microscopic theory, the fluid would tend to a normal 
one. 
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Table 1. The values of w and 7 corresponding to the best fit for a given choice of D. The error bars correspond to the 
projections of the 95% C.L. ellipses on the two axes. Also shown are the corresponding values of qo and z tr . 



As the x 2 -values listed in table 1 show, the fits are excellent. To further compare the shape of the 
theoretical spectrum with the Union2 data set, we also performed a Kolmogorv-Smirnov test. For each of 
the cases the K-S statistic was found to be smaller than 1.8 x 10~ 3 reflecting an extremely good fit. In fact, 
so good are the fits, that the current data is unable to differentiate between these choices of parameters. 
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Figure 3. 95% C.L. contours in the 7 — w plane for different values of D. The points represent the best fit values for 
each D. 

Since D assumes only discrete values, we refrain from treating it as a free parameter for the rest of 
the analysis. Rather, for a given D, we consider the 7-10 plane as a two-dimensional parameter space. We 
may, then, attempt to define 95% C.L. contours in this plane by considering A% 2 - These are displayed in 
figure 3. As is evident, there is a strong negative correlation between the two parameters. Note, furthermore, 
that positive w (hence, positive pressure) is essentially ruled out. Similarly, integral values of 7 are also 
essentially ruled out. Further, this is true even if we consider values of D far larger than those preferred by 
microscopic theories of high energy physics. As can be deduced from table 1 , each of the marked points in 
figure 3 denotes essentially a global minimum of x 2 , with the position of minima getting increasingly closer 
as one increases D arbitrarily. 
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Once a(r) and, hence, H(z) is determined in a model, one may also calculate both the present decel- 
eration parameter qo as well as zt r , the redshift corresponding to the epoch of transition from the decelerated 
to the accelerated phase. Note that these values are not uniquely determined by the data alone as the cosmo- 
logical model has a strong bearing on this determination. Also shown, in table 1, are the values of qo and 
as determined within our model. 
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Figure 4. Evolution of the deceleration parameter q with the redshift. The numbers in the parentheses refer to 

5 Conclusions and Discussion 

Several approaches have been explored in the literature to arrive at a late-time acceleration of the Universe. 
However, none of the models are completely satisfactory. Most models rely on a mysterious type of mat- 
ter (Dark Energy) which gives an effective repulsive gravity which is supposed to provide the observed 
accelerated expansion. This mysterious source term is perceived only through its effect on cosmological 
expansion. 

In this paper, we have followed an approach using higher dimensions. The mysterious behaviour of 
matter (a la the Dark Energy) is manifested directly only in the extra dimensions. The accelerated expansion 
of the scale factor in normal dimensions is produced only indirectly by the source term through its effect on 
the extra dimensions. Thus, in our model, the issue of invoking strange properties for the dark energy fluid 
does not arise as it exhibits its unusual property only in the extra dimensions. 

The particularly simple scenario we consider here is able to produce the requisite late time acceleration. 
The matter content acts as a pressureless gas in the normal dimensions and has a monomial equation of state 
as far as the extra dimensions are concerned. Just this simple ansatz leads to not only a late time acceleration, 
but also to a veiy moderate contraction of the extra dimensions since the epoch of radiation-matter equality 
although the normal dimensions have expanded by a factor of ~ 3000 in the same time period. With the 
size of the extra dimensions hardly changing, the role of any fields confined to the extra-dimensions as also 
those of any possible Kaluza-Klein towers of the SM fields in the low energy limit has remained essentially 
unaltered in this epoch. 

For a significantly wide range of parameters, the model shows excellent agreement with the obser- 
vational data on Type la supernovae, as attested to by both a low \ 2 P er degree of freedom as well as the 
Kolmogorov-Smirnov statistic. We still need a negative pressure nonetheless, albeit limited to the extra 
dimensions. Whether a more complicated scenario, involving a non-isotropic extra dimensional subspace 
and/or multiple fluids, obviates this restriction is yet to be seen. 

An interesting feature of our model is that the present phase of accelerated expansion is, generically, a 
transient one. In figure 4, we plot q(z) against log 10 (l+^). Note that the future is defined by log 10 (l + z) < 
0. For example, (D, w, 7) = (2, —3.71, 0.72), namely one of our best fit points, the universe will transit to 
a decelerating phase at z = —0.23. Similarly, for (D, w, 7) = (6, —2.80, 0.74), which is somewhat away 
from a best fit (but within the 95% C.L region), this transition would occur at z = —0.28. 

One aspect which we have not dwelt with in this paper is the rate of growth of perturbations. While 
a detailed analysis is beyond the scope of this paper and would be addressed elsewhere, let us make a few 
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comments. With the dynamics presented here being the dominant driving mechanism post matter-radiation 
equality, changes in structure formation would be confined to this epoch. In fact, the quantum of difference 
from the A-CDM scenario would be of the same approximate size as in a large class of theories with a 
dynamical Dark Energy source. This is supported by the fact that qo in this model (see table 1) agrees 
(within error bars) with —0.6, the value it assumes in the standard A-CDM scenario. Indeed, a careful 
analysis of this aspect could prove of value in further narrowing down of the parameter space, with larger D 
being less preferred. 
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